 ssgsea_score <- GSVA::gsva(
     expr = exp_normalized[, group_infor$sample] %>% as.matrix(),
     gset.idx.list = pathway,
     method = "ssgsea",
     ssgsea.norm = TRUE,
     verbose = T,
     parallel.sz = 50
 ) %>% as.data.frame() # signature 'matrix,list'